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ABSTRACT 

Jet systems with two unequal components interact with their parent accretion 
disks through the asymmetric removal of linear momentum from the star-disk 
system. We show that as a result of this interaction, the disk's state of least 
energy is not made up of orbits that lie in a plane containing the star's equator 
as in a disk without a jet. The disk's profile has the shape of a sombrero curved 
in the direction of acceleration. For this novel state of minimum energy, we 
derive the temperature profile of thin disks. The flaring geometry caused by the 
sombrero profile increases the disk temperature especially in its outer regions. 
The jet-induced acceleration disturbs the vertical equilibrium of the disk leading 
to mass loss in the form of a secondary wind emanating from the upper face 
of the disk. Jet time variability causes the disk to radially expand or contract 
depending on whether the induced acceleration increases or decreases. Jet time 
variability also excites vertical motion and eccentric distortions in the disk and 
affects the sombrero profile's curvature. These perturbations lead to the heating 
of the disk through its viscous stresses as it tries to settle into the varying state 
of minimum energy. The jet-disk interaction studied here will help estimate the 
duration of the jet episode in star-disk systems and may explain the origin of the 
recently observed one-sided molecular outflow of the HH 30 disk-jet system. 

Subject headings: accretion, accretion disks — stars: mass loss — stars: winds, 
outflows 



Introduction 



Stellar jets interact dynamically with their parent accretion disk in tw o fundamental 
ways . First, they extract angular momentum from the inner part of the disk (IPudritz et al. 
20061 ). a theoretical prediction that ha s been confirmed recently by the obse rvation of velocity 
asymmetries in the jet of RW Auriga (IWoitas et al.ll2005l ; iRay et al.ll2006l ). Second, a stellar 
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jet has been proposed as the physical carrier of hi gh temperature minerals from the Sun's 



vicinity to the o uter parts of the early solar system (jSkinnerlll990l ; iLiffman and Brownlll995 



Shu et al.lll996l ). The discovery of Calcium- Aluminum rich inclusions in the grains of Comet 
Wild 2 collected by Stardust seem to confirm this hypothesis. 

In this paper, we study a new dynamical interaction between the disk and jet that arises 
from the asymmetric removal of linear momentum from the star-disk system. A growing num- 



ber of stellar jets appear to be asymmetric as t 



differ by about a factor of 2 (IHirth et al. 



1994 



le ejection velocities of the jet and counterjet 



Woitas et al.ll2002l : lLopez-Martm et aUl2003 



Lavalley et al. 1997 ; Fernandez and Comeron 20051 ). Jet launching regions are confined to 



the inner part of the disk with estimates from 0.0 1 AU for the X-win d model to a few AU 
for disk-wind models (e.g. 1.6 AU for RW Auriga. I Woitas et al.ll2005l ). The integrated mo- 
mentum loss over the launching region accelerates the center of mass of the star-disk system 
which coincides with the star's center for axisymmetric systems. A gas element in the outer 
disk sees the matter around which it revolves accelerating and responds dynamically to such 
change. This response is the object of this work. 

An order of magnitude for the acceleration is obtained from the observed mass loss rates 
and ejection velocities in disk-jet systems as: 



A ~ 10 
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where M is the star's mass augmented by that of the je t launching region of the disk, M 
is the mass loss rate, and v e is the jet's ejection velocity (IPodio et al.ll2006l ). This estimate 
is an instantaneous lower bound on the total momentum loss as we lack long time span 
observations of stellar jets as well as observations of the jet engine within a few AU from 
the star. Although small, the acceleration amplitude (CQ) was shown to offer a possible 
explanation for the large eccentricities and secular res onances o f extra solar planets' orbits 
provided the planets formed in a jet-sustaining disk (jNamounil (120051 ) hereafter Paper I). 
Over the duration of the asymmetric jet episode, the star-disk system acquires a residual 
velocity, V, that must be smaller than the stellar velocity dispersion in the Galaxy, (v g ). 
This constraint yields an upper bound on the duration of the asymmetric jet episode, r, as: 



^ in5 3 x 10- 12 kms 2 (v g ) 

r < 10 5 N 91 years, 

A 10 kms" 1 



(2) 



where we used V ~ At. Asymmetric momentum removal would contribute at most a few 
kms" 1 to the observed ve locity dispersion of young stars (of order a few tens of kms -1 , 
Binney fc Merrifieldl (119981 )) leading to an upper bound on the jet's lifetime comparable to 
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the viscous time of the disk. The observation of (the same) asymmetric jets since 1994, 
yields a lower bound of ~ 10 years on how long asymmetric jets may be sustained. 

Jet asymmetry may have two origins: the a mbient medium and /or an asymmetric 
magnetic field that threads the disk. The former (IFerreira et al.l 120061 ) has the advantage 
of not adding more complex geometries to the magnetic field responsible for jet generation. 
The latter option of an asymmetric magnetic field has been considered in the context of AGN 
systems where it has been shown to warp the inner part s of the accretion disks through the 
asymmetric local magnetic pressure (jWang et al.lll992l ). Future observations inside a few 
AU from the central star will help elucidate the origin of jet asymmetry. 

In this work, we show that the acceleration of the star-inner-disk system with respect 
to the outer disk (section 2) alters the disk's state of minimum energy from a plane to a 
sombrero-shaped surface curved along the direction of acceleration (section 3). This flaring 
geometry has new consequences for the stellar irradiation of the disk (section 4) as well as for 
its vertical structure (section 5). Jet time- variability is reflected in the varying acceleration 
imparted locally to the disk leading to the global radial expansion and contraction of the 
disk (section 6). These results are summed up in section 7 by illustrating the disk's response 
as the jet is launched, and discussing their possible application to the HH30 disk-jet system. 



2. Jet-induced acceleration 



Jet-induced acceleration arises from the removal of momentum from the star-disk system 
(Fig. [TJ in an inertial frame. To see this, we divide the star-disk system into three particles: 
the star, an inner disk materialized by a particle that looses mass and linear momentum 
located at 0.5 AU, and the outer disk materialized by a constant mass particle located far 
outside, for instance at 50 AU. This is justified because the jet launching region is confined 
to the inner part of the disk. Gauss's theorem (the divergence theorem) for the gravitational 
potential of the form r -1 implies that the outer particle (outer disk) orbits around the 
combined masses of the star and the inner particle (the inner disk) which form what we 
shall call the central object. The outer particle is not concerned by the details of mass and 
momentum losses in the central object around which it revolves, it only responds to them. 
The central object or more accurately, the center of mass of the star and the inner disk, 
is loosing mass and linear momentum. It is therefore being accelerated with respect to the 
outer particle by an amount equal to the ratio of the momentum loss rate from the inner 
particle divided by the total mass of the central object -note that the inner particle itself 
is accelerated by the ratio of the momentum loss rate to its mass. The two-body problem, 
central object-outer particle, being symmetric, this also amounts to saying that the outer 
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particle is accelerated with respect to the central object that remains stationary. This shows 
that the outer particle evolution may be described by the equation: 

dv GM 

at |x| d 

where x and v are the position and velocity and G is the gravitational constant. The 
mass M is that of the star augmented by that of the inner particle and corresponds to the 
total mass contained in the jet launching region. The acceleration A is obtained as the 
ratio of the total momentum loss to the total mass M. In practice, we can neglect the 
contributions of the disk's mass and its variations in M for two reasons. First, the mass of 
the inner disk is small compared to the mass of the star. Second, the mass loss is small and 
amounts to a correspondingly small dynam ical effect g iven by the Jeans radial migration 



rate r/r = 10~ 8 yr _1 for M = 10~ 8 MqJt' 1 (IJeansI 119241 ) . 



From this simple picture, we can understand what goes on in an extended continuous 
disk. Mass and linear momentum loss modify the disk's equations of motion by including 
the disk's gravitational potential as well as the momentum loss from the disk (in the case 
of disk winds). The contribution of the disk's potential is negligible because the disk's mass 
is small and the migration resulting from mass loss is also small. As for the momentum 
loss, it is a function of the radius, r, that represents the integrated momentum loss up to r. 
Calling tjlr the outer radius of the jet launching region, then for r > tjlr, the acceleration 
is not a function of r as it represents the total momentum loss from the star-disk system. 
For r < tjlr, the variation with r of the momentum loss depends on the model of jet 
generation. Being interested in the flaring of the disk which mostly concerns its outer parts, 
we may employ Eq. (J3J) to describe the evolution of a gas element located outside of the jet 
launching region where A is independent of r. 



3. State of minimum energy: The sombrero profile 

We consider the steady state of the disk and accordingly take A to be constant. The 
dynamical problem described by equation d3J) has two constants of motion, the energy E = 
v 2 /2 — GM/\x\ — A • x, and the projected component of the angular momentum h = x x v 
on the direction of acceleration h z = h • A/A 

Before discussing the situation where the jet is not perpendicular to the disk, we first 
assume that the acceleration lies along the disk's plane normal. Least energy orbits make 
up the invariable surface of the disk - if there is no acceleration, this surface is the disk's 
midplane. We determine the orbits of least energy by making use of cylindrical coordinates 
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9, z) and by expressing the conservation of the vertical component of angular momentum 
= r 2 9. The energy equation then reads: 



E 



r 2 + z 2 



+ 



h 2 
2r 2 



GM 



Az. 



7 (4) 

Vr z + z z 

Minimizing the energy integral yields f = z = which corresponds to circular orbits and 
shows that r and z are not independent as the star's pull balances the acceleration A ac- 
cording to: 

GM7. 

A. (5) 



(j-2 _|_ z 2\3/2 

The vertical location of a circular orbit of radius r and its orbital rotation rate, Q 
given as: 



2 


GM 


, Q = 


(r 2 + z 2 f 2 



9, are 



(6) 



We see that orbits of least energy are circular but do not lie in the equator plane of the star, 
instead they hover above it in the direction of acceleration. The locus of circular orbits in 
the rz-plane is shown in Fig. ([2J upper panel). However, not all circular orbits are stable 
in the long term. There is a natural truncation radius around the star where the star's 
gravitational acceleration becomes comparable to A (Paper I). Beyond this radius, the star's 
pull is weak and orbits escape its gravity. The location of this radius denoted ctk p ir is given 
by the equality of the acceleration excitation time Ta = 27rrSl/3|v4| (Paper I) and the orbital 
period T = 2ir/Q. It is given as: 
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(7) 



An accretion disk accelerated by a jet uses the stable orbits of the minimum energy state. 
The profile thus obtained is that of a sombrero as shown in Fig. lower panel). We call 
z = S(r) the stable disk profile obtained by inverting Equation ([6]). We will also use the 
linearized profile, S(r) for z r which reads: z = Ar 3 /GM = T 3 /3a£ plr . We remark that 
the sombrero profile coincides with the steady state profile of a thin jet-sustaining disk. For 
(realistic) thick disks, the sombrero corresponds to the invariable surface, the equivalent of 
the disk's midplane if there were no acceleration, as well as to the disk's surface of maximum 
density (section 5). In the case where the acceleration is not initially perpendicular to 
the disk, the state of least energy is still a sombrero whose symmetry axis is that of the 
acceleration. To reach this state, the disk experiences a transient phase during which it 
dissipates energy (section 6). 
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4. Temperature 



The stellar irradiation of a flat thin disk (i.e. a planar disk with zero thick ness) with 
no je t-induced acceleration yields a decreasing temperature profile with distance (IFriedjung 
19851 ). It is given as T d oc T # (R*/r) 3 / 4 where T d is the disk's temperature, T* and R* are the 
star's temperature and radius respectively. For a non-planar disk, the disk's temperature at 
distance r from the star is related to the stellar temperature through: 



T 4 



rp4 



rz' - z 2R* 



2r 



3nr 



R* 



-i 2 



(IKenyon and Hartmannl 119871 ) where z is the vertical location of the disk and z' = dz/dr. 
For z < r, we substitute the linearized sombrero profile z = S(r) into (jHJ) to find: 



Ta — T* 



Rl 



3a kplr 



+ 



2^ 

3nr 2 



(9) 



The outer regions (r ^> R*) are those that contribute most to the temperature profile. 
Accordingly, the disk's temperature is given as: 



T H ~ 13 



4000 K 



2R Q 



500 AU 

flkplr 



K. 



(10) 



The fundamental state of the disk therefore has a residual temperature regardless of radius. 
The sombrero profile cancels out the r~ 2 distance factor responsible for the decrease of 
temperature with distance. Relaxing z r and using the full shape of S(r) leads to a 
modest temperature increase with distance. The increase and decrease of the slopes of z/r 
on both sides of the sombrero for a finite thickness disk will make the temperature grow 
beyond the estimate of Eq. ffTOl for the upper side of the disk and exposes the lower side to 
some irradiation. The unequal heating of jet-sustaining disks therefore sets a local vertical 
temperature gradient that leads to the stratification of the disk. 



5. Hydrostatic equilibrium and wind generation 

For a disk with a finite thickness, we may seek the vertical profile obtained from hydro- 
static equilibrium. The hydrostatic equilibrium equations are now modified by the presence 
of the vertical acceleration A. Assuming an isothermal disk of pressure p = c 2 p where p is 
the density and c the sound speed, the hydrostatic equilibrium equation reads: 

dp GMzp 



dz ( r 2 + ^2)3/2 



P A = 0, (11) 
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whose solution is 



P 



Po{r Jexp — 



Vr 2 + z' 



+ Az 



(12) 



The sombrero profile here corresponds to the surface of maximum density. With respect 
to this profile and where z -C r, the disk's thickness, H = c/fl, is similar to that of an 
unaccelerated disk. This solution however has a flaw: the linear term related to the accel- 
eration yields an infinite disk mass when p is integrated over all space. The finite value of 
p for large positive z is indicative of a matter outflow from the upper side of the disk in 
response to the lowering of the vertical pull of the star by the acceleration A. Consequently, 
Eq. (TTTT) can no longer describe the vertical state of the disk. The steady state properties 
of th e new wind component m ay be described by the unmagnetized Grad-Shafranov equa- 
tion (IScott and Lovelace! 119821 ) subject to the finite duration acceleration A. We expect the 
density and velocity profile of the wind component to bear some resemb lance to the hydro- 
dyna mic jets generated as force-free vortex funnels in accretion disks (jScott and Lovelace 



1983). 



6. Variability 



6.1. Vertical evolution 



S tellar jets are time - variable processes (IRaga fc Kofmanl Il992l ; lHartigan fc Raymond 



19931 : iFalle fc Ragal Il995l : iReipurth et al.l 120021 : iRav et al.1 12006^ To ascertain the conse- 
quences of the jet-induced acceleration's time variations, we first consider an initially circular 
planar disk subject to a vertical acceleration of magnitude A. This models an instantaneous 
increase of acceleration from zero to A. Solving the vertical equation of motion in ([3]) under 
the assumption that z <Cr yields: 



Ar 3 
GM 



cos Qt). 



(13) 



The vertical motion is an oscillation with respect to the linearized sombrero profile S(r) 
on a timescale equal to the local orbital period in the disk, T. The oscillation amplitude is 
given by the distance from the initial state which in this case is the equator plane to the 
location of the sombrero profile. If the acceleration grew from zero to A on a timescale, 
|yl/v4|, longer than the orbital period, T, the disk's midsurface height z would follow the 
evolving sombrero profile adiabatically and there will be no oscillation with respect to the 
surface of least energy. To illustrate the two excitation regimes, sudden and adiabatic, 
we show in Fig. ([3]) the response of the disk's midsurface height obtained by numerically 
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integrating the equation of motion ([3]). The acceleration's time dependence is taken as 
A = An ax (l — exp[— t/100yr]) x W(t) where A max = 8 x 10~ 12 kms~ 2 , the same magnitude 
as in Fig. (j2J), and W(t) is a window function that vanishes outside the range [0,5000 yr] 
and equals unity inside it. The minimum truncation radius is therefore 500 AU. As the 
launching timescale is = 100 yr, the disk evolves adiabatically to the sombrero profile 

inside 22 AU as can be seen on z(t)/r at r = 10 AU (T = 32 yr) in Fig. Outside 
22 AU, the dynamical time T is longer than the launching timescale of 100 yr implying that 
jet launching is perceived as a sudden excitation by most of the outer disk whose typical 
response is described by Eq. (|13p and shown in Fig (j3J) for r = 100 AU. 



The acceleration is shut off suddenly after 5000 yr. This determines an outer limit, 
r ~ 300 AU, in the disk beyond which the state of least energy remains the disk's initial 
state (star's equator plane) as the dynamical time outside this radius is too long to travel all 
the way to the location of the sombrero profile (Eq. [131) . Beyond this radius, the jet-induced 
acceleration episode is felt as a short (with respect to the orbital period) velocity pulse to 
the outer disk resulting in the excitation of vertical motion with respect to the disk's initial 
midplane as seen in Fig. (jHJ) for r = 400 AU. This type of excitation has been discussed 
extensively in Paper I (section 6). The sudden absence of acceleration is felt everywhere 
in the disk, where for instance at r = 10 AU, the lack of the jet-induced acceleration that 
balanced the star's vertical gravitational pull leads to a vertical oscillation with respect to 
the equator plane and with an amplitude equal to the distance from the sombrero to that 
plane. This effect is the same as the excitation described by Eq. (|13|) only this time the 
invariable surface is the disk's midplane and the initial state is the sombrero profile. For 
both the accelerated and unaccelerated phases, the vertical oscillations with respect to the 
surface of least energy (be it a sombrero if A ^ or a plane if A = 0) are damped by 
the viscous stresses of the disk. For instance, in a more realistic simulation taking into 
account the disk's viscosity, the observed oscillations about the sombrero height z/r ~ 10~ 2 
for r = 100 AU, would be erased and the disk would settle into the sombrero profile before 
the acceleration is shut off at t = 5000 yr. We note that the illustration of the disk's 
midsurface height response to the acceleration's growth shown here is not peculiar to the 
launching of the jet and equally represents the disk's response to more general variations of 
the jet-induced acceleration as those associa ted with a varying ejection velocity at the source 



(IRaga fc Kofmanlll992l ; iFalle fc Ragalll995( ). 
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6.2. Horizontal evolution 

For a time-variable acceleration, the conservation of the specific angular momentum 
h z = r 2 Q forces the mean orbital radius to change in order to compensate for the increase of 
the midsurface height z to the location of sombrero profile. The rate of this change is found 
using the conservation of h z only as: 

r = 3 - z, (14) 
r 

where we neglected a term of second order in z/r. The rate (fl~4"l) is valid regardless of the 
time dependence of A and whether the initial state of the disk is planar. The change in 
radius is of second order in A and leads to a radial expansion of the disk if the acceleration 
increases with time. If however A decreases, the midsurface height z is lowered to the 
new sombrero location and Eq. ( Till) shows that the disk undergoes a radial contraction. 
The expansion and contraction occur on the local dynamical time T if it is larger than the 
variability timescale (as in Eq. [TBI) . Otherwise, radial motion occurs adiabatically 

with respect to variability. Outside the radius where T = it is the former evolution 

that takes place implying that contraction and expansion start from the inside of the disk 
and proceed outward (as T increases with r). The local radial displacement Ar is given as: 

|Ar| _ 3eA 2 r 4 _ e 
~V ~ G 2 M 2 ~ 3 

where e = 1 if the acceleration's change is sudden, < T. For an adiabatic change, 

|yl/v4| ^> T, the averaging of the linearized equations of motion ([3J yields e = 2/3. The 
outer half of the disk is the most affected by the expansion and contraction: at the disk's 
median radius, r = ak p i r /2, \Ar\/r ~ 2%. We also note that this rate is much larger than 
that of the Jeans migration (section 2). 

Figure (jlj) shows the local expansion and contraction of the disk at different radii under 
the same conditions as those of the previous section and Fig. ([3]). At r = 10AU the expan- 
sion is adiabatic as the launching timescale |^4/y4| is larger than the dynamical time, while 
the contraction is sudden as the acceleration is abruptly turned off. Where the dynamical 
timescale is much larger than the acceleration's duration as for r = 400 AU, there is a net 
radial expansion outward as the whole acceleration episode is felt as a velocity kick to the 
outer disk (Paper I, section 6) — note for instance how for r = 400 AU, the expansion 
appears to be delayed until the end of the acceleration phase. 

Another consequence of the acceleration's variability is the local excitation of eccentricity 
in the disk. Using the linearized equations of motion ([3]) for the radius r, it can be shown 



a kplr 



(15) 
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that the local eccentricity growth rate due to the vertical acceleration is: 

3A 2 r 4 Q 



71 



G 2 M 2 ' 



(16) 



with a secular turnaround time equal to The maximum eccentricity is therefore e max = 
( r /ok P ir) 2 /3. We note that the eccentricity gradient d(re)/dr = (r/a^) 2 is always smaller 
than unity implying that the forced eccentric distortions of the gas streamlines do not make 
them cross and consequently they do not shock. The estimate of Eq. ( Fl6|) is valid during 
the acceleration phase. Outside the location where the dynamical time is smaller than the 
duration of the acceleration phase, the net effect of acceleration amounts to a small amplitude 
velocity kick imparted to the disk. To the fluid elements of the disk's midplane, the velocity 
kick imparts an eccentricity given as (Paper I, section 6): 

V 2 r/GM 

eW = l + l4/GM - (17) 

where V is the total velocity obtained by integrating the acceleration over its duration. The 
global response of the disk's viscous stresses to the eccentric disturbance (fl6j) depend s on the 



presc ription of the viscous stresses that is employed as well as on the density profile ( jOgilvie 



200ll ). We expect the viscous stresses to damp the eccentric distortions as well as the vertical 
oscillations leading to the heating of the disk. In the outermost part of the accretion disk 
where the dynamical time is long compared to the acceleration's duration and where the 
viscous timescale is similarly long, we expect jet-sustaining accretion disks to exhibit global 
eccentric distortions. 



7. Synthesis 

Stellar jets are ubiquitous processes of star formation and star-disk interaction. The 
asymmetric removal of linear momentum from the star-disk system accelerates the outer 
part of the disk with respect the jet launching region. The acceleration amplitude is small 
([T]) and can be modeled as a time dependent vector directed along the star's rotation axis. 
The state of minimum energy of the disk is a sombrero-shaped surface made up of circular 
planar orbits that hover above the star's equator plane ([6]) increasingly higher with distance. 
As the jet is launched, the initially circular equator disk moves towards the state of minimum 
energy and the truncation radius moves from infinity to the location given by dTJ). If the 
disk's size is larger than the minimum truncation radius over the acceleration phase, the 
disk will be truncated by the jet-induced acceleration. There is a specific radius in the disk 
where the orbital period T matches the jet launching timescale Inside this radius, the 
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disk moves away from the equator plane to the instantaneous sombrero surface adiabatically. 
Outside the matching radius, vertical motion ffTB"]) as well as eccentricity distortions (fT6j) with 
respect to the sombrero are excited. The conservation of angular momentum implies that 
as the jet is launched the disk expands radially (I14|ll5j) . Viscous stresses damp the vertical 
motion and the eccentricity distortions to heat up the disk especially in its outer regions. 
The steady state of the disk thus assumes a sombrero profile that exposes the upper face of 
the disk to the star's light canceling at first order the irradiation's decrease with distance 
f llOp . The disk cannot maintain hydrostatic equilibrium under the jet acceleration (TT2"j) . 
Instead the disk looses mass from its upper side giving rise to a secondary asymmetric wind 
component. When the jet-induced acceleration decreases, the sombrero looses curvature, the 
disk contracts radially at the rate given by (fT5l) and the truncation radius moves outward 

The mechanisms discussed above only require the onset of asymmetric momentum re- 
moval in the star-disk system. Whether momentum is removed from the star itself through 
some pec uliar stellar wind com ponent, or more realistically from the disk or the star-disk 
interface (jFerreira et al.l 120061 ) . the outer disk will respond to the integrated momentum 
loss from the matter around which it gravitates. The mechanisms discussed above also do 
not have a time requirement on the acceleration's duration. The evolution towards the som- 
brero, the radial contraction and expansion proceed from the inner part of the disk outwards. 
Therefore, these mechanisms affect the disk inside the radius where the dynamical time is 
comparable to the acceleration's duration. In this sense, disk observations can be used to 
constrain the duration of the jet episode and/or asymmetric momentum removal by iden- 
tifying the outermost radius where relaxation toward the sombrero profile has taken place. 
Asymmetric momentum removal imprints two more signatures on the star-disk system: the 
presence of a global eccentric distortion in the outermost part of the disk, and an enhance- 
ment of the random component of the stellar velocity in the Galaxy. The eccentric distortion 
results from a perceived sudden pull of the star-inner disk system with respect to the outer 
disk along the jet system's axis and should be an increasing function of radius (fTTj) that does 
not depend on the acceleration's timescales. The change in the stellar velocity dispersion 
may be ascertained by statistically analyzing young stars' proper motion at different epochs 
of their evolution. 



A good illustration of some of the new propertie s discussed above is found in the HH 
30 disk-jet system. This system ( jBurrows et al.lll996l ) is a nearly edge-on 420 AU disk that 
appears as two reflection nebulosities separated by a dark lane. The jet and counterjet 
appear to be asymmetric although both the ejection speed and the mass loss rate remain 
uncertain (jBurrows et al.l Il996l ; iBacciotti et al.lll999f). The j et exhibits variability in the 
form of knots with periods of 2.5 years (jBurrows et al.lll996l ). The disk's variable lateral 
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asymmetries appear mainly on one side of the disk with possible periods fro m a few days to 
a few years Jstapelfeldt et alill999l : IWatson et aPl200ol : ICotera et allbpOlh . More recently, 
Pety et al.l ( 20061 ) reported the observation of a molecular outflow emanating from the upper 
nebulosity. The outflow was traced back to a scale smaller than the size of the disk while 
no rotation signatures were detected above 1 kms -1 at 200 AU. The variable lateral asym- 
metries are currently explained by the presence of a single hot spot on the star requiring 
a com plex magnetic field stru cture that differentiates substantially between the two stellar 
poles (jStapelfeldt et al.l Il999l ) . If the one-sided outflow resulted from the entrainment of 
surrounding material then the cloud structure would be substantially different on both sides 
of the disk. In the context of the jet-disk interaction studied here, it is natural to attribute 
the lack of lateral asymmetries in the lower nebulosity to a change in the disk's mean pro- 
file, and to identify the one-sided molecular outflow as the new secondary wind component. 
Further modeling of the disk temperature, wind density and velocity profiles produced by 
the jet-induced acceleration will help determine the properties of the HH 30 disk and its 
molecular outflow. 



The author thanks the referee for constructive comments that helped improve the clarity 
of the paper and Christiane Froeschle for stimulating discussions and the careful reading of 
the manuscript. This work was supported by the Programme National de Planetologie. 
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Fig. 1. — Schematic of the star-disk system. The momentum loss from the upper component 
is larger than that of the lower component. The net momentum loss from the jet launching 
region, which extends to a few astronomical units from the star, modifies the dynamics of 
the outer disk according to Eq. 
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Fig. 2. — The locus of circular orbits in the rz-plane (AU) for an acceleration A = 8 x 
1CT 12 kms~ 2 . The truncation radius a kp i r = 500 AU. The upper panel shows all possible 
orbits: stable (solid line) and unstable (dashed line). The lower panel is a zoom of the stable 
orbits inside a radius of 200 AU around the star (dashed circle not to scale). 
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Fig. 3. — Vertical evolution of the outer disk. The equations of motion ([3]) were integrated 
numerically to simulate the vertical motion of a gas element initially located in the disk's 
midplane before the jet-induced acceleration is activated. The acceleration is modeled as 
A = 8 x 1(T 12 kms~ 2 x (1 - exp[-t/100yr]) x W(t) where W(t) = 1 if < t < 5000 yr and 
zero otherwise. The four panels (left scale) illustrate z(t)/r for r =10 AU, 100 AU, 250 AU 
and 400 AU -the corresponding dynamical times T are given in each panel. The dashed-line 
profile (right scale) in the four panels shows the acceleration A normalized to its maximum 
value 8 x 10~ 12 kms~ 2 . 
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Fig. 4. — Radial expansion and contraction of the outer disk. The equations of motion ([3]) 
were integrated numerically to simulate the radial motion of a gas element initially located 
in the disk's midplane before the jet-induced acceleration is activated. The acceleration 
is modeled as A = 8 x 1(T 12 kms" 2 x (1 - exp[-t/100 yr]) x W(t) where W(t) = 1 if 
< t < 5000 yr and zero otherwise. The four panels (left scale) illustrate Ar/r for r =10 AU, 
100 AU, 250 AU and 400 AU -the corresponding dynamical times T are given in each panel. 
The dashed-line profile (right scale) in the four panels shows the acceleration A normalized 
to its maximum value 8 x 10 -12 kms~ 2 . 



